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Effects  of  the  Hubbard  interaction  and  electrostatic  pinning  in  polyacetylene 
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The  discrete  model  for  the  coupled  electron-phonon  system  in  polyacetylene  has  been  extended  to 
include  the  electronic  effects  of  electrostatic  pinning  and  the  Hubbard  interaction  in  the  unrestrict¬ 
ed  Hartree-Fock  approximation.  To  fit  experimental  structural  data,  the  Hamiltonian  also  includes 
the  lattice  kinetic  energy  and  the  coupling  of  several  vibrational  degrees  of  freedom.  This  model  is 
then  used  to  self-consistently  determine  the  electronic  eigenstates  and  the  structural  configuration 
for  an  odd-membered  ring.  The  electronic  states  from  this  equilibrium  configuration  are  then  used 
to  calculate  the  charge  and  spin  density,  optical  absorption,  and  the  single-electron  loop  corrections 
to  the  dynamical  matrix.  The  resulting  vibrational  modes  are  examined  and  used  in  conjunction  . 
with  the  self-consistent  electronic  states  to  calculate  the  infrared-absorption  spectra.  Assuming  that 
the  pinning  arises  from  an  oppositely  charged  soliton  on  an  adjacent  chain  leads  to  a  pinning  fre¬ 
quency  of  =5(X)  cm  ',  in  excellent  agreement  with  observed  results.  A  value  for  the  Hubbard  in¬ 
teraction  of  4  eV  leads  to  an  excellent  fit  for  the  rest  of  the  infrared  spectra  and  the  optical  absorp¬ 
tion  for  both  a  charged  and  a  neutral  soliton. 


1.  INTRODUCTION 

The  simple  model  of  polyacetylene,  (CH)^,  originally 
proposed  by  Su,  Schrieffer,  and  Heeger'  (SSH),  has  been 
very  successful  in  explaining  many  experimental  results. 
In  order  to  develop  a  more  quantitatively  accurate  model 
it  is  important  to  consider  additional  interactions  outside 
the  range  of  the  original  model.  Thus  in  this  study  we 
propose  to  investigate  the  effects  of  both  the  Hubbard  in¬ 
teraction  and  electrostatic  pinning  on  the  electronic  and 
optical  (visible  and  infrared)  properties  of  polyacetylene 
in  the  presence  of  a  soliton.  The  effects  of  pinning  from  a 
point  charge  on  the  optical  absorption  in  (CH)j,  has  been 
considered  for  a  variety  of  configurations  by  Glick  and 
Bryant.^  The  “minimum  energy  configuration"  in  that 
calculation  was  not  obtained  from  a  minimum  energy 
self-consistency  relation,  nonetheless  they  obtained  a 
rather  nice  fit  with  experimental  results  for  the  optical 
absorption  of  a  charged  soliton  well  removed  from  the 
end  of  the  chain.  In  this  work  we  consider  the  pinning 
resulting  from  a  point  charge,  a  reasonable  point  of  view 
for  some  chemically  doped  samples,  and  in  an  effort  to 
understand  the  finite  "pinned”  frequency  for  optically  ex¬ 
cited  (CH),,  the  pinning  resulting  from  an  equal  and  op¬ 
positely  charged  soliton  on  an  adjacent  chain.  (We  have 
recently  become  aware  of  two  independent  studies  using 
this  same  pinning  concept,  one  using  a  phenomenological 
approach,^  and  another  which  calculates  numerically  the 
phonon  eigenstates  in  the  single  component  problem  *  to 
predict  the  frequency  of  the  “pinned"  mode  in  optically 
excited  samples.) 

The  effects  of  the  Hubbard  interaction  in  the  Hartree- 
Fock  approximation  on  the  electronic  states  in  polyace¬ 
tylene  when  localized  defects  are  present  have  been  con¬ 
sidered  with  some  success  both  numerically  in  the 
discrete  model’'*  and  analytically  in  the  continuum  mod¬ 
el.^  With  this  in  mind  we  have  decided  to  incorporate  the 
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same  approximation  to  calculate  the  electronic  states  in 
the  discrete  model  for  polyacetylene.  We  then  extend 
these  calculations  to  obtain  the  optical  absorption  for 
both  charged  and  neutral  solitons,  the  phonon  modes  for 
charged  solitons,  and  the  infrared  absorption  in  polyace¬ 
tylene.  Recently  the  phonon  modes  with  particular  in¬ 
terest  in  the  bound  states  have  been  studied  in  the  pres¬ 
ence  of  the  Hubbard  term*  and  a  combination  of  the 
Hubbard  term  and  pinning.^  The  effects  of  pinning  and 
electron-electron  repulsion  appear  to  do  a  reasonable  job 
of  fitting  the  frequencies  of  the  “pinned”  mode  for  both 
chemically  doped  and  optically  excited  samples.^  How¬ 
ever,  since  the  vibrational  modes  in  polyacetylene  are 
more  complex  than  those  in  a  simple  one  dimensional  sin¬ 
gle  component  model  it  is  necessary  to  extend  this  model 
to  one  which  includes  three  structural  degrees  of  free¬ 
dom.  A  very  useful  scheme,  developed  by  Horovitz’ 
(amplitude-mode  or  AM  formalism)  has  been  used  with 
the  framework  of  the  continuum  model'®  to  study  the 
spectrum  of  modes  with  large  oscillator  strengths  associ¬ 
ated  with  the  translation  of  mobile  defects,  “7"  modes. 
However  it  has  been  shown  that  additional  ir  active 
bound  modes  in  the  single  component  problem  can  result 
in  additional  satellite  features  (/4  ~)  in  the  infrared  spec¬ 
tra."''^  Thus  the  existence  of  the  ir  active  third  bound 
mode  for  a  single  structural  degree  of  freedom  in  both  the 
continuum  model”  and  the  discrete  model'*" for  po¬ 
lyacetylene  should  result  in  additional  features  in  a  realis¬ 
tic  calculation  of  the  infrared  absorption.  In  fact  these 
additional  features  have  been  predicted  in  a  multicom¬ 
ponent  calculation  in  the  absence  of  pinning  and  Hub- 
herj  interactions  in  both  the  continuum  model'*  and  the 
discrete  model  with  a  semicontinuum  approximation.'* 
Calculating  the  phonon  modes  and  the  ir  absorption  in 
such  a  multiple  component  model  enables  one  to  not  only 
fit  the  experimental  frequencies  and  oscillator  strengths. 
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of  the  three  infrared  active  “7"  modes  which  have  rela¬ 
tively  large  large  oscillator  strengths  but  the  frequencies 
and  oscillator  strengths  of  the  satellite  A  localized 
modes  as  well.  The  prediction  of  these  additional  ir  ac¬ 
tive  modes  has  stimulated  some  recent  experimental  in¬ 
terest  in  the  weak  satellite  (A  infrared  absorption 
peaks.^'*’^'  The  measured  oscillator  strengths  for  these 
A  ^  modes  is  smaller  than  that  predicted  by  coupling  the 
phonon  modes  to  the  electron-phonon  induced  current 
fluctuations  in  both  the  continuum  model  calculations'* 
and  the  discrete  model  calculations.''*  We  will  show  that 
with  pinning  and  a  value  of  4  eV  for  the  Hubbard  interac¬ 
tion  an  excellent  fit  with  infrared  absorption  data  is  ob¬ 
tained  for  both  (CH)_j  and  (CD)_,. 

Additionally  recent  work  by  the  authors  has  indicated 
that  the  width  of  the  pinned  mode  may  be  an  inherent 
property  of  polyacetylene  due  to  strong  mixing  of  the  low 
frequency  translational  mode  with  low  frequency  acoustic 
phonons  creating  hybrid  states  in  the  acoustic  phonon 
branch  with  large  infrared  oscillator  strengths. This 
would  eliminate  the  need  for  a  distribution  of  pinning  po¬ 
tentials  in  the  material.  However,  we  show  in  this  study 
that  as  the  soliton  becomes  pinned  this  mixing  is  dramati¬ 
cally  reduced  and  as  originally  postulated  a  distribution 
of  pinning  potentials  or  pinning  parameters'*  is  required 
_ I 

//  =  —  ^  +  \  ,a^n,a  ^n,<T^n  +  l.a  ^  "T  S  +  I 

n,(j  ^  n 


to  explain  the  experimentally  observed  width  in  the 
“pinned  mode." 

11.  ELECTRONIC  EIGENSTATES 

To  model  polyacetylene  we  assume  the  Su,  Schrieffer, 
and  Heeger  (SSH)  Hamiltonian*  with  on  site  repulsion  of 
the  Hubbard  type  and  screened  Coulomb  electrostatic 
pinning.  The  electrostatic  pinning  is  only  considered  for 
charged  systems  and  is  modeled  as  either  a  point  charge 
off  the  chain  or  as  a  distributed  charge  from  an  equal  and 
opposite  soliton  residing  on  an  adjacent  chain.  To  model 
the  vibrational  modes  we  extend  the  SSH  Hamiltonian  to 
include  multiple  structural  degrees  of  freedom  in  a  matter 
analogous  to  that  developed  for  the  continuum  model  by 
Horovitz’  (AM  formalism).  In  all  cases  considered  we  in¬ 
itially  calculate  a  self-consistent  structural  configuration 
and  the  corresponding  eigenstates  which  depend  strongly 
upon  the  electrostatic  pinning  and  the  strength  of  the 
Hubbard  term.  Since  we  are  primarily  concerned  with 
the  effects  on  the  localized  electronic  state  that  occurs  in 
this  problem,  we  will  consider  the  Hubbard  term  in  the 
mean  field  unrestricted  Hartree-Fock  approximation 
when  calculating  the  eltotronic  eigenstates. 

The  SSH  Hamiltonian  with  Coulombic  pinning  and  on 
site  Hubbard  repulsion  is  given  by 


2 


(2.1a) 


with  the  hopping  energy  being  proportional  to  the  spac¬ 
ing  between  the  carbon  sites, 

‘n.n.,l  =  h-Y(yn+\-yJ-  <2.  lb) 

Tne  c’s  create  and  annihilate  ir  electrons  on  the  appropri¬ 
ate  site  and  y  is  the  electron-phonon  coupling  constant. 
U  is  the  Hubbard  constant  with  V(n)  being  the  site 
dependent  Coulomb  pinning  potential,  and  K  represents 
the  a  electron  spring  constant  of  the  polyacetylene  chain. 
In  the  unrestricted  Hartree-Fock  approximation  the  on 
site  Hubbard  repulsion  is  replaced  with 

U'L[P-o(n)-\]UlaCn.a-l^-  f2.1c) 

n,a 

The  quantity  p_Jn)-\  is  simply  the  excess  charge  with 
negative  spin  where  is  given  by 

<2.2) 

k 

and  2*  is  sum  only  over  occupied  states.  We  choose 
K(/j)  to  be  a  screened  electrostatic  Coulomb  potential  in 
an  anisotropic  dielectric  media  of  the  form 

,  -Fe'  ^  [pJm)-\]exp(~AR„,„) 

F(«)= - 7=-  2 - ^ - .  <2.3a) 

(iV  m.o  ^n,m 

where 


R 


2 

n.  m 


in  ^  ^ 

ft 


(2.3b) 


For  pinning  in  optically  excited  samples,  we  assume  from 
x-ray  diffraction  data  for  polyacetylene^^  that  a,  the  dis¬ 
tance  between  carbon  sites^  is  1.23  A  and  d,  the  inter¬ 
chain  separation  is  4.24  A.  Numerical  values  of  the 
dielectric  constant  both  parallel  and  perpendicular  to  the 
chain  have  been  determined  experimentally^'*  to  be 
fi  =  7. 1  and  ^2=1- 8-  The  integer  m  sums  over  the  car¬ 
bon  sites  on  an  adjacent  chain  and  pjm)  —  j  is  taken  to 
be  the  negative  of  the  chain  under  consideration.  T  is  a 
parameter  that  chosen  to  be  2  to  take  into  account  the  re¬ 
duced  mass  effect  for  a  pair  of  solitons.  A  was  chosen  to 
fit  the  experimental  data  and  found  to  be  0.59.  For  pin¬ 
ning  in  chemically  doped  systems  the  expression 
2„[p„(m)—  |]  is  replaced  by  a  Kronecker  delta  inunction 
which  centers  a  point  charge  at  m  =(fi  -|- 1  )/2,  and  F  for 
this  case  is  1.  Since  polyacetylene  is  known  to  expand 
when  it  is  doped, the  distance  of  the  charged  impurity 
from  the  chain  is  taken  to  be  0.52d  rather  than  0.5d.  In 
Sec.  V,  where  we  discuss  the  infrared  activity  of  polyace¬ 
tylene,  it  will  be  seen  that  this  choice  of  parameters  does 
an  excellent  job  of  fitting  the  frequencies  of  the  "pinned” 
modes  in  polyacetylene. 

The  above  Hamiltonian,  Eq.  (2.1),  yields  in  the  unre¬ 
stricted  Hartree-Fock  approximation  the  following 
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difference  equation  for  the  electronic  eigenstates, 

^k.a'l'kj  l.n'PkJ  "  -  1  >  +  rt  +  1  ) 

+  \  (2.4) 

Thus  the  eigenstates,  and  eigenvalues,  „,  are  ob¬ 
tained  by  simply  numerically  diagonalizing  this  expres¬ 
sion.  To  consider  the  problem  of  a  chain  with  a  single 
soliton,  the  total  number  of  sites,  N,  is  always  chosen  to 
be  odd,  and  to  avoid  difficulties  witli  end  effects  periodic 
boundary  conditions  are  assumed.  Taking  a  variational 
derivative  of  this  Hamiltonian  with  respect  toy^  leads  to 
the  self-consistency  relation, 

k,a 

-f  -f  1  )-t- const.  (2.5) 

To  model  polyacetylene  we  have  chosen  2/o  =  5  eV, 
‘*)'3'o/2tQ  =  Ao/2fo=0. 14  (  2Ao  being  equal  to  the  Peierls 
gap)  which  leads  to  a  value  of  the  dimensionless 
electron-phonon  coupling  constant  of  K  =  Sy^/rrK{2lQ) 
=  0.4197. 

This  self-consistency  expression  combined  with  Eq. 
(2.2)  for  both  possible  spins  enables  all  of  the  input  pa¬ 
rameters  to  Eq.  (2.3)  to  be  determined  from  its  own  solu¬ 
tions.  Thus  Eq.  (2.3)  can  be  diagonalized  repeatedly  with 
new  values  of  y„  +  |  — y,  and  pjn}  calculated  for  each  di- 
agonalization  until  a  minimum  energy  self-consistent 
solution  is  obtained  within  a  specified  convergence  cri¬ 
teria.  From  the  resulting  electronic  eigenstates  it  is 
straightforward  to  determine  charge  and  spin  densities 
for  charged  and  neutral  solitons  as  well  as  such  parame¬ 
ters  as  the  soliton  excitation  energy,  gap  state  energy  lev¬ 
el,  and  the  soliton  width  as  a  function  of  the  Hubbard  U. 
Quantitative  values  for  these  parameters  can  be  very  in¬ 
formative,  for  example  the  ratio  of  spin  densities  on  even 
and  odd  sites  {p../p^)  for  a  neutral  soliton  is  strongly 
dependent  on  the  strength  of  the  Hubbard  term. 
Electron-nuclear  double-resonance  (ENDOR)  measure¬ 
ments  have  determined  this  ratio  to  be  — From  our 
calculations  this  is  consistent  with  U  =  3  —  4  eV,  which 
is  consistent  with  earlier  results.’"^ 

The  full  Hubbard  term  is  given  by 

^Huh  =  Y  j).  (2.6) 

i,c 

If  we  define  n,  =n,  and  the  site  dependent  spin 

density  s,  =  jin,  ^  —  then  we  can  rewrite  the  Hub¬ 

bard  term  as 


"Hub  =  t/2 

I 


(n,-  l)(n,  -n 

4 


(2.7) 


For  a  charged  soliton  is  zero,  thus  to  minimize  the  en¬ 
ergy  the  soliton  width  increases  while  for  a  neutral  soli¬ 
ton  H;  —  I  is  zero,  and  the  effect  of  the  Hubbard  term  is  to 
decrease  the  soliton  width.  This  is  seen  in  our  results 
where  the  relative  soliton  width  increases  from  1  to  1.18 


for  a  charged  soliton  and  decreases  from  1  to  0.82  for  a 
neutral  soliton  as  U  ranges  from  0  to  4  eV.  Again  this  is 
consistent  with  results  of  other  investigators.’"’  Also 
consistent  with  this  interaction  is  a  shift  in  the  gap  state 
energy.  For  a  charged  soliton  we  find  that  the  gap  state 
energy  shifts  from  =0  for  U  =0  to  £^^^  =  ±0. 134  eV 
for  t/  =4  eV  with  the  -f  iign  corresponding  to  an  n-type 
soliton  and  the  —  sign  corresponding  to  a  p-type  soliton. 
For  a  neutral  soliton  the  energy  shifts  of  the  midgap  state 
are  more  significant  and  the  spin  degenv.iai.>  is  lificd.  As 
1/  ranges  from  0  to  4  eV  the  gap  state  for  spin  up  (the  oc¬ 
cupied  state)  shifts  from  Eg^p  =  0  to  £gjp=— 0.396  eV 
with  the  corresponding  shift  up  for  the  spi  down  gap 
state.  This  will  have  some  significance  in  the  optical  ab¬ 
sorption,  as  we  will  see  in  the  next  section. 

With  the  solutions  to  the  self-consistent  Hamiltonian 
described  above,  we  are  able  to  calculate  many  interesting 
properties  of  polyacetylene.  In  addition  to  the  charge 
and  spin  density  discussed  above,  we  use  these  states  to 
find  the  optical  absorption,  single  electron  loop  correc¬ 
tions  to  the  dynamical  matrix  and  the  corresponding  vi¬ 
brational  states,  and  the  infrared  absorption  spectra. 
Each  of  these  calculations  will  be  discussed  in  the  follow- 
ing  sections. 


in.  OPTICAL  ABSORPTION 

Since  we  are  dealing  with  an  extended  system,  it  is 
awkward  to  deal  directly  with  induced  moments;  conse¬ 
quently,  we  consider  current-current  correlations  to  cal¬ 
culate  the  conductivity,  o(<u).  With  the  Hamiltonian  of 
Eq.  (2.1),  the  current  density  is  expressed  as 

J(r)  =  t  +  —  (3.1) 

n,(7 


In  the  Matsubara  formalism  the  current -current  correla¬ 
tion  is  given  by 

rr{T)=-{T,J(T)J(0)),  (3.2) 

where  is  the  imaginary  time-temperature  ordering 
operator.  The  retarded  current-current  correlation  func¬ 
tion  is  obtained  by  Fourier-series  expanding  in  Matsu¬ 
bara  frequencies  and  then  analytically  continuing 
— ►w-fib.  The  frequency  dependent  conductivity  is 
derived  from  this  as 


(7(w)=  Re 


i  ir(  fc)  -b  iS ) 
0) 


(3.3> 


Following  the  usual  formalism,  where  Wick’s  theorem  is 
used  to  expand  the  current-current  correlation  function 
in  terms  of  Green’s  functions  and  then  using  an  eigen¬ 
function  expansion  for  these  Green’s  functions,  yields  in 
the  low-temperature  limit 


<^(«)  =  -  T' fk.,. of,. k.a  6(tu^ (3.4a) 


(O 


k.q.a 


with 


fk.q.o  =  ^  +n0^,<,(l)-tl';k,„(l)tfr,,,,(/  -1-1)], 


(3.4b) 
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where  sums  over  empty  k  states  and  filled  q  states. 

To  be  entirely  consistent  with  our  approximation  of  the 
effect  of  the  Hubbard  term  on  the  optical  absorption,  we 
should  include  additional  effects  to  the  current-current 
correlation  other  than  direct  transitions  between  renor- 
mali2ed  states.  The  result  of  such  an  effort  would  result 
in  summing  over  a  frequency  dependent  matrix  element 
times  the  same  joint  density  of  states  delta  function.  This 
would  be  only  a  small  correction  to  optical  conductivity 
and  would  significantly  increase  the  difficulty  of  the  cal¬ 
culation  and  hence  is  ignored. 

The  results  of  the  calculations  are  somewhat  dependent 
on  the  width  of  the  Lorentzian  used  to  represent  the  delta 
function  in  sum  above.  We  typically  use  as  small  a  value 
as  possible  and  still  have  the  absorption  appear  as  a 
smooth  curve.  For  a  101  site  chain  this  requires  the  full 
width  at  half  maximum  (FWHM)  to  be  0.06(2fo)-  The 
effect  on  the  optical  absorption  of  the  Hubbard  term  for 
charged  solitons  is  not  surprising  in  that  the  main  effect 
is  simply  a  small  shift  in  the  midgap  absorption  band. 
However,  as  mentioned  above  the  shift  is  more  dramatic 
for  a  neutral  soliton,  and  for  the  same  Lorentzian  half- 
width  the  inner-band  transitions  merge  with  interband 
transitions  at  about  U  =4.0  eV.  In  Fig.  1  we  plot  the  op¬ 
tical  absorption  in  a  neutral  soliton  for  both  {/=3  eV 
and  U  =A  eV.  Since  the  experimentally  observed  optical 
absorption  spectra  are  broader  than  those  plotted  in  Fig. 
1,  we  would  expect  from  these  calculations  that  the  inner 
band  transition  could  not  be  resolved  experimentally  for 
a  value  of  U  between  3  eV  and  4  eV.  Thus  a  value  of  U 
between  3  eV  and  4  eV  is  again  consistent  with  experi¬ 
mental  observation,  where  very  careful  measurements  of 
optical  absorption  in  ammonia-compensated  polyace¬ 
tylene  show  no  sign  of  inner  band  transitions.^^ 


FIG.  1.  Optical  absorption  for  a  neutral  soliton  with  (dashed) 
f/=3.0eVand  (solid)  (/=4.0eV. 


IV.  VIBRATIONAL  MODES 


To  make  direct  experimental  comparison  of  structural 
vibrations  we  have  chosen  to  extend  the  multicomponent 
concept  of  the  AM  formalism’  to  the  discrete  SSH  model 
for  polyacetylene.  This  model  Hamiltonian  includes  the 
lattice  kinetic  energy  and  couples  the  n  electrons  to  a  set 
of  displacement  fields,  as  well  as  the  Hubbard  and 
pinning  terms, 


^  2  "  -Hi  ^  -r  \  ,a^n,a  "h  +  l.a  ^  “h  2  2  +  1  3'a.n  ^  +  2  2  ^  a. 1 

a,n  a,n 


(4.1) 


To  maintain  a  fixed  electronic  gap,  in  analogy  with  the 
multicomponent  continuum  model,  we  take 
2oya3'o,n=m.  thus  leaving  ,  unchanged.  From 
self-consistency  we  find  (^a/Ta'Ta.n  iTn*  and 

defining  a  dimensionless  electron  phonon  coupling  for 
each  component  as  kg  =  %y\/trKa(2tQ)  leads  to 
'LaK=^- 

In  a  system  with  strong  electron-phonon  coupling  such 
as  this,  it  is  necessary  to  consider  the  screening  induced 
by  the  electron-phonon  interaction  via  charge  fluctua¬ 
tions  in  the  tt  electron  gas.  This  is  conveniently  done  us¬ 
ing  single  electron  loop  perturbation  theory  depicted  di- 
agrammatically  in  Fig.  2(a).  The  adiabatic  solution  of  the 
integral  equation  represented  by  this  diagram  in  a  mul- 


€>  =  O  ^  0--0 


(C)  - - - -t- -  0-- 

FIG.  2.  The  diagrammatic  representations  of  (a)  the  integral 
equation  for  the  vibrational  modes,  (b)  electron-phonon  induced 
charge  fluctuations  including  all  first  order  contributions  from 
the  Hubbard  interaction,  and  (c)  the  integral  equation  for  the 
screened  Hubbard  interaction. 
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ticomponent  representation  is  equivalent  to 

,  <"a)' 

o)lk<l>k^ct,n)=—^[2<i>,_t,{a.n)-4>,  i,(a,n  +  1  *(a,n  -  D] 

o  A.  A 

p,m  I 


(4.2) 


where  IKn.m)  is  the  electron-phonon  induced  polariza¬ 
tion  represented  by  the  shaded  bubble  in  Fig.  2(a).  In  this 
expression  (<u°  )^  is  the  bare  k  =0  optical  phonon  frequen¬ 
cy  squared  and  is  given  by  and  the  correspond¬ 

ing  dimensionless  coupling  constant  A.„  =  8y3/7r/fg(2ro). 
These  constants  have  been  determined  uniquely  from  Ra¬ 
man  scattering  data^*  and  are  listed  in  Table  I  for  both 
(CH)^  and  (CD),.  In  a  multicomponent  description 
both  the  coupling  constants,  A.„,  and  bare  frequi  ncies,  <u®, 
will  be  k  dependent.  However,  the  vibrational  modes  of 
interest  are  all  near  or  at  k  =0  on  either  the  acoustic  or 
optical  phonon  branch.  Since  the  electron-phonon  in¬ 
teraction  leaves  the  vibrational  modes  in  the  acoustic 
phonon  branch  virtually  unchanged,  and  since  the  pa¬ 
rameters  are  taken  from  a  continuum  model  fit,  which  is 
appropriate  for  k  =0  optical  phonons,  we  assume  that  ig¬ 
noring  the  k  dependence  of  these  constants  will  still  result 
in  physically  meaningful  results. 

Since  we  are  considering  a  many  body  term  (Hubbard 
interaction),  the  outer  electron  lines  in  the  polarization 
bubble  are  bold,  indicating  that  they  include  the  first  or¬ 
der  contribution  of  the  Hubbard  interaction  to  the  proper 
self-energy.  In  a  many  body  formalism  this  usually  re¬ 
quires  the  solution  of  a  Dyson  equation;  however,  for  this 
particular  problem  that  is  identical  to  the  solution  of  the 
unrestricted  Hartree-Fock  problem  as  described  in  Sec. 
II.  To  be  consistent  with  that  level  of  approximation  in 
our  calculation  of  11,  we  consider  the  additional  contribu¬ 
tions  of  the  Hubbard  interaction  to  the  electron-phonon 
induced  polarization  shown  in  Fig.  2(b),  where 
II  =  ll'"-)-n'^'.  As  shown  in  Figs.  2(b)  and  2(c),  this  ad¬ 
ditional  contribution  requires  an  infinite  sum  of  bubble 
diagrams  connected  with  the  Hubbard  constant  U.  Thus 
to  solve  the  integral  equation  (4.2)  above  first  requires 
summing  all  the  contributions  to  the  polarization.  We 
perform  this  sum  by  first  solving  the  integral  equation 
that  is  represented  diagrammatically  in  Fig.  2(c).  If  we 


I 

define  in  the  adiabatic  approximation  the  bubble  with 
two  Hubbard  vertices  [Fig.  2(c)]  as  n'®’(n,n'),  and  the 
dressed  Hubbard  interaction  as  then  the  integral 

equation  in  Fig.  2(c)  is  written  as 

yn.n=^^n.n-+^  1  Il'°V  «,  m  )  (4.3a) 


where 


Il'®'(n,n')  =  22 

k.g 


(4.3b) 


Since  we  intend  to  use  these  states  to  calculate  the  in¬ 
frared  absorption  in  polyacetylene,  we  are  only  interested 
in  the  phonon  modes  for  a  charged  soliton.  For  the  case 
of  a  charged  soliton  the  electronic  eigenstates  are  spin  de¬ 
generate,  thus  we  have  dropped  their  spin  dependence 
and  included  a  factor  of  2  in  Eq.  (4.3b).  To  solve  this  in¬ 
tegral  equation  for  we  start  by  simply  diagonalizing 
n'®'(M,/j').  Il‘®’  can  then  be  expanded  in  terms  of  its 
orthonormal  set  of  eigenfunctions  as 

(4.4) 

k 


In  terms  of  these  states  the  solution  for  V„  „■  is  simply 


k 


<f>k{rt)<f>k{n') 


(4.5) 


The  second  diagram  in  Fig.  2(b),  which  we  define  as 
although  somewhat  tedious,  is  now  straight¬ 
forward  to  evaluate.  Keeping  in  mind  that  the  outer  ver¬ 
tices  involve  matrix  elements  of  the  electron-phonon  in¬ 
teraction  and  that  the  internal  vertices  involve  matrix  ele¬ 
ments  of  the  Hubbard  interaction,  this  diagram  in  the 
adiabatic  approximation  is  expressed  as 


n'^’(n,n’)  =  4  2  2 

m.m'  k,q 


Ek-E,- 


(4.6) 


TABLE  I.  Coupling  constants  and  bare  optical  phonon  fre- 
quencies  for  (CH),  and  (CD),- _ 


'I 

— 

£i>"  (cm 

') 

a 

(CH), 

(CD). 

(CH), 

(CD). 

1 

0.07 

0.06 

1234 

921 

2 

0.02 

0.005 

1304 

1207 

3 

0.91 

0.935 

2040 

2040 

where  the  electron-phonon  vertex,  g*,,(n),  is  given  by 

gt.,(n)  =  V'*(")[<A,(«  -  \)] 

-f- t/',(/l)[  !/’*(«  -I- 1  )-!/>*(«  —  1 )].  (4.7) 

Again  with  the  factor  of  4  we  have  taken  advantage  of 
the  spin  degeneracy  of  the  electronic  eigenstates.  The  ex¬ 
pression  for  the  first  diagram  on  the  right-hand  side  of 
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FIG.  3.  The  optical  (upper  curves)  and  acoustic  (lower  curves) 
components  of  the  translation  phonon  mode  for  (dashed)  U  =0 
eV  and  (solid)  U  =4  eV. 


the  equals  sign  in  Fig.  2(b)  is  the  usual  one  in  the  SSH 
model  for  adiabatic  electron-phonon  screening  of  the  vi¬ 
brational  modes, 


k.q 


(4.8) 


[["'  and  11'^'  are  then  summed  and  Eq.  (4.2)  is  diagonal¬ 
ized  to  determine  the  phonon  eigenstates  of  the  mul¬ 
ticomponent  model. 

As  expected  from  continuum  model  results,  this  exten¬ 
sion  of  the  discrete  model  to  include  several  structural 
degrees  of  freedom  results  in  three  ir  active  localized 
modes  related  to  the  translational  mode  of  the  soli- 
tof,s.  11,12  three  additional  ir  active  modes  related  to 
the  third  bound  mode  in  the  single  component  version  of 
the  SSH  model. In  the  absence  of  pinning  the  effects 
of  the  Hubbard  term  have  relatively  little  effect  on  these 
bound  modes.  However,  the  coupling  of  the  translational 
mode  to  the  low  frequency  acoustic  modes  is  enhanced  as 
C/  increases  from  0  eV  to  4  eV.  This  is  apparent  in  Fig.  3, 
where  for  simplicity  we  plot  the  acoustic  and  optical 
components  of  the  translational  mode  for  the  single  com¬ 
ponent  model  for  both  (/  =0  eV  and  2/  =4  eV.  The.se  re¬ 
sults  are  physically  reasonable,  as  the  Hubbard  term 
causes  a  local  expansion  of  the  lattice  in  the  region  of  the 
.soliton.  The  mixing  for  both  cases  is  significantly  larger 
than  the  original  calculation,'^  which  used  a  semicontinu- 
um  model^'*  solutions  for  the  electronic  eigenstates.  Thus 
if  this  mixing  were  maintained  in  the  presence  of  pinning, 
it  could  explain  the  width  of  the  pinned  mode  as  an  in¬ 
herent  property  of  the  system. 

The  inclusion  of  pinning  has  a  rather  dramatic  effect 
on  the  frequency  of  the  translational  mode  as  well  as 
some  smaller  shifts  in  the  higher  frequency  modes.  With 
the  choice  of  pinning  parameters  discussed  and  in  the  ab¬ 
sence  of  the  Hubbard  term,  the  frequency  of  the  pinned 
mode  for  an  optically  excited  sample  is  shifted  from  near 
zero  to  530  cm  '  and  for  a  chemically  doped  .sample  to 
805  cm^'  for  (CH), .  For  (CD)^  the  corresponding  fre¬ 


quency  shifts  are  to  493  cm  ’  and  778  cm  ‘,  respective¬ 
ly.  These  frequencies  are  somewhat  high  (low)  for  the  op¬ 
tically  excited  (chemically  doped)  polyacetylene.  Howev¬ 
er,  when  we  include  a  Hubbard  interaction  with  U  =4 
e\',  these  frequencies  are  shifted  for  both  the  optically 
and  chemically  doped  systems.  For  (CH)^  and  (CD)^  the 
pinned  frequencies  for  optically  excited  samples  become 
503  cm''  and  474  cm"',  respectively,  while  for  chemi¬ 
cally  doped  systems  the  corresponding  frequencies  are 
877  cm“'  and  767  cm' '.  These  results  are  in  excellent 
agreement  with  experimental  results,  implying  that  the 
mechanisms  proposed  in  Sec.  II  for  pinning  are  physical¬ 
ly  reasonable.  However,  at  the  pinned  frequencies  the 
wave  number  of  the  acoustic  phonons  is  much  larger 
than  the  inverse  of  the  characteristic  width  of  the  pho¬ 
nons,  thus  the  mixing  of  the  translational  mode  with 
these  states  is  dramatically  reduced.  This  reduces  the 
transfer  of  ir  oscillator  strength  to  the  nearby  acoustic 
modes,  creating  a  relatively  sharp  ir  active  mode  at  the 
pinned  frequency.  Consequently,  the  experimentally  ob¬ 
served  width  of  the  pinned  mode  is  almost  certainly  due 
to  a  distribution  in  pinning  potentials  (or  pinning  param¬ 
eters)  in  the  polymer,  and  not  the  result  of  a  distribution 
of  oscillator  strengths  due  to  mixing  with  acoustic  pho¬ 
nons. 


V.  INFRARED  ABSORPTION 

In  Sec.  Ill  we  discussed  the  optical  absorption  from  the 
point  of  view  of  retarded  current-current  correlations. 
This  bare  contribution  (with  respect  to  the  electron- 
phonon  interaction)  to  the  optical  activity  is  represented 
by  the  single  loop  in  Fig.  4(a)  which  contributes  only  the 
allowed  electronic  transitions  to  the  optical  conductivity. 
Turning  on  the  lattice  kinetic  energy  couples  the  ionic 
motions  to  the  charge  fluctuations  of  the  a-  electron  gas. 
This  coupling  transfers  oscillator  strength  from  the  pure¬ 
ly  electronic  transitions  to  lower  frequency  modes  associ¬ 
ated  with  the  IR  active  vibrational  modes.  The  relevant 
correction  is  expressed  diagrammatically  in  Fig.  4(b), 


(b) 


'  O '  o  ^  c>—o 

FIG.  4.  The  diagrammatic  representations  of  (a)  the  current- 
current  correlation,  (b)  correction  to  current-current  correlation 
that  includes  low  frequency  ir  active  modes,  .and  (c)  current- 
current  correlation  including  all  first  order  contributions  from 
the  Hubbard  interaction. 
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which  makes  a  contribution  to  the  current-current  corre-  ming  all  relevant  bubble  diagrams  connected  by  the  Hub- 
lation  function  that  oscillates  at  the  frequency  Con-  bard  term  as  shown  in  Fig.  4(c)  and  use  the  adiabatic  ap- 

sistent  with  our  previous  efforts  in  calculating  the  vibra-  proximation  where  appropriate. 

tional  modes,  we  self-con.sistently  consider  the  effects  of  The  real  part  of  the  infrared  conductivity  for  several 
the  Hubbard  term  on  the  current  fluctuations  by  sum-  lattice  degrees  of  freedom  is  expressed  as 


<7(«)=Re  —  2  +  D„p{m,m‘;ct)  +  ib)  y^{m',o)  +  iS) 

^  m. a,0 


where  the  Green’s  function,  D^  p,  representing  the  vibra¬ 
tional  modes,  is  given  by 

<(>, 

(co  +  ib}^--ct)^^ 

(5.2) 

The  coupling  between  the  vibrational  modes  and  the  elec¬ 
tron  gas  is  mediated  by  the  electron-phonon  interaction 
whose  effects  on  the  dynamical  problem  were  discussed  in 
the  last  section.  The  lowest  order  contribution  (with 
respect  to  the  Hubbard  interaction)  to  the  current  fluc¬ 
tuations,  j°{m,a)  +  ib),  is  linear  in  co  and  given  by 


g{m,m'-,co  +  ib)=—  - -  T 

VM^Mp  ft 


/( m, w -H i6 )  =  2"  fk.qSg.k^  fn ) 


2(6; -I- <5 ) 

(o}  +  ib)^-(E^-E.)^ 


where  for  the  spin  degenerate  system  being  considered, 
^  is  simply  twice  the  original  definition  of  f  k,g,a-  The 
higher  order  contributions  to  j  are  found  with  a  tech¬ 
nique  exactly  analogous  to  that  described  in  Sec.  IV.  As 
an  example  we  note  that  the  expression  representing  the 
last  diagram  in  Fig.  4(c)  is  nearly  identical  to  Eq.  (4.6); 
the  only  change  is  ^  replacing  the  first  gt,,(n).  The 
effective  coupling  constant,  y„/M^,  is  more  conveniently 
expressed  in  terms  of  the  dimensionless  coupling  constant 
an^the  bare  phonon  frequency  in  the  functional  form 
l/ A.„6;°.  Since  the  vibrational  frequencies  are  small  com¬ 
pared  to  the  electronic  gap,  the  adiabatic  approximation 
for  this  case  is  equivalent  to  ignoring  only  the  frequency 
dependence  in  the  denominator,  and  in  general  the  ex¬ 
pression  for  j  is  linear  in  o)  and  of  the  form 
j \  m,<o)  =  2(oK'(m).  The  expression  for  the  low  frequency 
optical  activity  now  reduces  to 


Thus  the  contribution  to  the  infrared  conductivity  is 
large  only  for  vibrational  eigenstates  which  overlap 
strongly  with  the  electron-phonon  induced  current  fluc¬ 
tuations.  In  general  these  fluctuations  are  localized  at  the 
defect  site  and  thus  only  localized  states  with  the  correct 
symmetry  can  exhibit  any  oscillator  strength. 

We  perform  the  dot  product  in  Eq.  (5.4)  numerically 
for  a  variety  of  Hubbard  interaction  strengths  for  both 
unpinned,  and  optically  and  chemically  pinned  solitons. 
The  results  for  both  (CH),  and  (CD)_,  are  very  enlighten¬ 
ing.  In  Fig.  5  we  present  the  results  for  unpinned  (CH)jt 
with  both  f/  =0  eV  and  C/  =4  eV.  The  most  noticeable 
change  is  the  increased  width  of  the  low  frequency  ab¬ 
sorption  due  to  mixing  with  the  acoustic  phonons. 
Another  significant  difference  is  the  increase  in  relative 
oscillator  strength  in  the  satellite  A  ~  modes.  However 
for  U  =0  eV  the  oscillator  strengths  in  these  modes  are 
substantially  larger  than  that  observed  experimentally 
and  of  even  more  significance  is  that  the  center  of  the 
pinned  mode,  depending  on  the  pinning  mechanism,  ex¬ 
ists  at  either  «500  cm“'  or  =i;850  cm~'.  To  better  ap¬ 
proximate  this  experimental  data  and  hopefully  under- 


tti(cm‘') 


FIG.  5.  The  infrared  absorption  in  (CH),  with  no  pinning  and 
a  Hubbard  term  of  (dashed)  {/  =0  eV  and  (solid)  t/  =4  eV. 
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stand  the  source  of  the  pinning,  in  Figs.  7(a)  and  7(b)  we 
plot  the  results  for  pinning  due  to  an  equal  and  opposite 
soliton  on  an  adjacent  chain  for  both  U  =0  eV  and  U  =4 
eV.  Other  than  the  shift  of  the  pinned  mode  to  a  finite 
frequency,  =530  cm“’,  the  satellite  peaks  have  virtually 
disappeared  for  f/  =0  eV.  Thus  pinning  depletes  the  os¬ 
cillator  strength  of  the  A  ~  modes  while  increasing  the 
oscillator  strength  of  the  modes  related  to  the  uniform 
translational  mode  of  the  soliton.  In  view  of  recent  stud¬ 
ies,  we  know  that  the  Hubbard  interaction  deepens  the 
effective  well  for  the  bound  state  vibrational  modes.*  The 
electrostatic  pinning  is  of  the  opposite  sign  of  the  Hub¬ 
bard  term,  making  it  reasonable  that  the  bound  state 
phonons  related  to  the  third  bound  mode  of  the  single 
component  problem  become  more  extended  and  carry 
less  oscillator  strength  (OS).  Additionally  the  transfer  of 
OS  to  the  translation  type  modes  is  consistent  with  the 
AM  formalism’  which  transfers  ir  activity  to  these  modes 
as  the  pinning  parameter  is  increased.  It  is  interesting, 
however,  that  a  value  f/  =4  eV  restores  a  significant 
amount  of  ir  activity  to  the  A  '  modes,  and  the  solid 
curve  in  Fig.  6(b)  does  an  excellent  Job  of  fitting  the  ex¬ 
perimental  data  for  optically  excited  (CH),,^*'^°  which  is 
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FIG.  6.  The  infrared  absorption  in  (CH),  with  optically  in¬ 
duced  pinning  for  (a)  U  =0  eV  and  (bl  U  =4  eV,  the  inset  here 
shows  experimental  data  for  the  frequency  range  1200  cm  '  to 
2000  cm  '. 


displayed  as  an  inset  in  this  same  figure.^’’  Note  the  sig¬ 
nature  of  the  A  “  in  both  the  experimental  and  theoreti¬ 
cal  results,  as  a  shoulder  on  the  dominant  absorption 
peak  at  1360  cm^  '.  When  including  pinning  and  a  (/  of  4 
eV,  we  obtain  frequencies  of  503  cm  ',  1280  cm*',  and 
1360  cm  “  '  for  the  translational  modes.  Their  relative  os¬ 
cillator  strengths  (OS)  are  44,  1.0,  and  5.0  respectively, 
which  compare  quite  favorably  with  measured  values.^' 
The  satellite  peaks  now  occur  at  in74  cm  ’  '  and  1456 
cm~'  with  OS  of  0.11  and  0.15.  For  completeness  and 
for  comparison  with  experimental  data.  Fig.  7  shows  the 
results  of  the  infrared  absorption  for  chemically  doped 
(CH),  with  U  =A  eV,  where  the  predicted  pinned  fre¬ 
quency  is  877  cm  '  with  an  OS  of  54.  The  additional 
two  translational-type  modes  occur  at  1283  cm  '  and 
1389  cm"'  with  relative  oscillator  strengths  of  1.0  and 
14.8.  At  first  glance  it  appears  the  OS  of  the  mode  at 
1283  cm"'  is  inconsistent  with  the  mechanism  of 
transferring  oscillator  strengths  to  these  modes  as  the 
pinning  is  increased.  This  can  be  easily  explained;  one  of 
the  A  "  modes  occurs  at  1280  cm"'  a.id  experimentally 
cannot  be  resolved  from  the  mode  at  1283  cm."  '  Since 
pinning  reduces  the  OS  of  the  A  "  modes,  the  net  effect  is 
a  slight  reduction  of  the  OS  of  this  peak.  In  Fig.  8,  we 
plot  the  results  for  (CD),  with  the  Hubbard  constant 
U  =4.0  eV  for  both  optically  and  chemically  generated 
solitons.  For  optically  excited  samples  the  pinned  mode 
occurs  at  474  cm"'.  The  additional  strongly  ir  active 
iranslational  modes  for  optically  excited  solitons  occur  at 
1042  and  1215  cm"'.  These  modes  have  OS  of  32.4,  9.8, 
and  1.0,  respectively,  which  again  are  in  relatively  good 
agreement  with  experimental  data.^'  For  these  condi¬ 
tions  the  only  observable  satellite  peak  occurs  at  1358 
cm"'  with  an  OS  of  0.15.  For  chemically  doped  (CD), 
we  predict  that  the  pinned  frequency  occurs  at  766  cm”' 
and  the  two  additional  translational  modes  at  1 138  cm  "  ' 
and  1233  cm"'.  These  modes  have  relative  oscillator 
strengths  of  2.0,  2.1.  and  1.0,  respectively.  Again  our  re¬ 
sults  predict  only  one  observable  satellite  peak,  which 
occurs  at  1402  cm "  '  with  an  OS  of  0.03. 

The  experimental  observation^'  of  the  satellite  feature 
in  (CH),  at  =  1050  cm  ~ '  is  somewhat  questionable  due 


FIG.  7.  Ttie  infrared  absorption  in  (CH),  with  chemically  in¬ 
duced  pinning  for  U  =  4  eV. 
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FIG.  8.  The  infrared  absorption  in  (CD),  with  (dashed)  opti¬ 
cally  and  (solid)  chemically  induced  pinning  for  t/  =4  eV. 


to  the  existence  of  a  large  naturally  occurring  ir  mode  at 
approximately  the  same  frequency  and  the  corresponding 
difficulty  of  accurately  subtracting  this  large  peak.  How¬ 
ever,  the  peak  at  approximately  1440  cm*'  has  no  such 
problem  and  is  a  clear  marker  for  A  "  vibrational  modes 
in  (CH),.  Since  pinning  is  required  to  fit  the  pinned 
mode  in  this  material,  then  as  we  have  seen  that  the  OS 
of  this  satellite  feature  at  1441  cm"'  can  be  used  to  quali¬ 
tatively  fit  the  strength  of  the  Hubbard  interaction  in 
(CH)_,.  Consistent  with  the  results  of  fitting  the  spin  den¬ 
sity  and  optical  absorption  for  a  neutral  soliton,  we  find 
that  a  good  fit  is  obtained  with  (f=;4  eV.  These  secon¬ 
dary  features  in  iCD)^  also  show  a  reasonable  fit  for  the 
same  value  of  U. 

VI.  CONCLUSION 

In  this  study  we  have  considered  the  effects  of  pinning 
and  the  Hubbard  interaction  on  various  electronic  and 
optical  properties,  as  well  as  an  extension  of  the  discrete 


SSH  model  to  include  several  degrees  of  freedom.  Several 
conclusions  are  possible  from  fits  with  experimental  data. 
First,  for  optically  excited  samples  our  postulated  mecha¬ 
nism  for  pinning,  an  equal  and  oppositely  charged  soliton 
on  an  adjacent  chain,  does  an  excellent  job  of  fitting  the 
frequency  of  the  pinned  mode  in  these  materials.  As  one 
would  expect  for  chemically  doped  samples,  a  point 
charge  off  the  chain  does  a  good  job  of  creating  a  poten¬ 
tial  well  with  the  correct  pinning  frequency,  providing 
that  the  distance  from  the  chain  is  increased  to  account 
for  some  expansion  upon  doping. We  have  also  demon¬ 
strated  that  the  width  of  the  experimentally  observed 
pinned  mode  is  probably  not  an  inherent  property  of  the 
system  due  to  a  transfer  of  CiCillator  strengths  to  modes 
in  the  acoustic  phonon  branch.  It  appears  that  a  distri¬ 
bution  of  pinning  potentials  (or  pinning  parameters)  is 
necessary  to  explain  this  width. 

Probably  the  most  interesting  result  of  this  study  is  the 
number  of  experiments:  spin  density  distribution  from 
ENDOR  measurements,^*  optical  absorption  for  a  neu¬ 
tral  soliton,^’  and  the  infrared  oscillator  strengths  of  the 
satellite  A  "  modes^' —  that  can  all  be  explained  by  in¬ 
cluding  an  effective  Hubbard  term  with  t/=:3  — 4  eV. 
The  consistency  of  all  of  these  results  validates  the  use  of 
the  unrestricted  Hartree-Fock  approximation  of  the  Hub¬ 
bard  term  when  considering  effects  of  electron-electron 
repulsion  for  a  chain  where  the  electronic  and  optical 
properties  are  dominated  by  the  localized  gap  state  that 
occurs  in  this  system. 

Comparisons  of  experimental  data  to  the  theoretical 
SSH  model  yields  excellent  agreement  when  we  addition¬ 
ally  include  (i)  electrostatic  pinning  for  charged  solitons, 
(ii)  the  Hubbard  term,  and  (iii)  several  structural  degrees 
of  freedom. 
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